Supersolid phase induced by correlated hopping in spin-1/2 frustrated quantum magnets 
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We show that correlated hopping of triplets, which is often the dominant source of kinetic energy in dimer- 
based frustrated quantum magnets, produces a remarkably strong tendency to form supersolid phases in a mag- 
netic field. These phases are characterized by simultaneous modulation and ordering of the longitudinal and 
transverse magnetization respectively. Using Quantum Monte Carlo and a semiclassical approach for an ef- 
fective hard-core boson model with nearest-neighbor repulsion on a square lattice, we prove in particular that 
a supersolid phase can exist even if the repulsion is not strong enough to stabilize an insulating phase at half- 
filling. Experimental implications for frustrated quantum antiferromagnets in a magnetic field at zero and finite 
temperature are discussed. 
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The identification of exotic states of quantum matter in mi- 
croscopic models is an important issue in current research on 
strongly correlated quantum systems. The existence of a su- 
persolid (SS) which simultaneously displays crystalline order 
(solid) and long-range phase coherence (superfluid, SF) has 
been definitely established recently thanks to extensive Quan- 
tum Monte Carlo (QMC) simulations of a hard-core boson 
model with nearest-neighbor (n.n.) repulsion on the triangu- 
lar lattice in the context of cold atoms loaded into optical lat- 
tices. [1, 2, 3] The possibility to realize a supersolid phase 
in dimer-based quantum magnets, first pointed out by Momoi 
and Totsuka in the context of SrCu2(BC<3)2 [4], has been fur- 
ther investigated very recently [5, 6, 7]. It relies on the de- 
scription of a polarized triplet on a dimer as a hard-core bo- 
son, a convenient language we will mostly use throughout this 
paper. Some trends have already emerged. In particular, for 
a supersolid to be realized, triplets induced by the magnetic 
field should have a small kinetic energy as compared to their 
mutual repulsion. In SU(2) models, this can be achieved if 
the inter-dimer coupling is frustrated, as noticed in the con- 
text of magnetization plateaux early on [8], and pointed out 
recently for supersolids by Sengupta and Batista [6]. This is 
not the whole story however. First of all, even if the kinetic 
energy is small, the transition between insulating and super- 
fluid phases can be first order, thus preempting the possibility 
of a supersolid phase. This is for instance the case of hard- 
core bosons with n.n. repulsion on the square lattice [9], in 
contrast to the results of Refs. [1, 2, 3] on the triangular lat- 
tice. Besides, when the hopping of a single triplet is strongly 
suppressed by frustration, the kinetic energy is not always sup- 
pressed accordingly. Indeed, correlated hopping terms much 
larger than single particle hopping are often generated. To see 
this, suppose for instance that dimers are locally arranged as 
in Fig. 1 with J' <C J, a situation typical of frustrated mag- 
nets. The amplitude for a triplet to hop from dimer 1 to 2 is 
strictly equal to 0, while the amplitude to hop from 1 to 3 is 
equal to J /2 /2J if dimer 2 is occupied by a triplet, as first 
pointed out in Ref. [10]. In two-dimensional arrangements 
of frustrated coupled dimers, single particle hopping does not 




t'=J'/2J 





FIG. 1 : Difference between single and correlated triplet hopping in 
a frustrated geometry, a) Spin language: thick solid (dashed) lines 
stand for dimer triplet (singlet), thin solid lines for inter-dimer cou- 
pling, b) Bosonic language: filled (open) circles denote hardcore 
boson sites which are occupied (empty). 



strictly vanish, but only appears at much higher order in per- 
turbation (6th order for the Shastry-Sutherland model realized 
in SrCu2(BC<3)2), while correlated hopping is typically a sec- 
ond order process [11]. So the problem of hard-core bosons 
with correlated hopping calls for further investigation. 

In this Letter, we show that supersolidity is strongly favored 
in dimer-based frustrated magnets because correlated hopping 
is the principal source of kinetic energy. The bottom line of 
our analysis can be summarized as follows: Under the effect 
of correlated hopping, a triplet cannot delocalize if it is alone, 
but it can in a crystalline arrangement of triplets with the ap- 
propriate geometry. In particular, if an insulating phase with 
a geometry compatible with correlated hopping is realized at 
commensurate filling, one may expect that upon adding par- 
ticles, the crystalline order will be retained, leading to a su- 
persolid phase. As we shall see, this simple mechanism gives 
rise to a tendency toward supersolidity which is much stronger 
than anticipated, in fact so strong that a supersolid phase can 
exist even if the corresponding solid phase cannot be stabi- 
lized. 

These conclusions are based on an extensive investigation 
with Stochastic Series Expansion (SSE) [12] Quantum Monte 
Carlo (QMC) simulations and with a semiclassical approxi- 
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FIG. 2: (Color online) Zero-temperature phase diagram for t = 0.95 
as a function of the 1/V versus fJ,/V. Open squares (closed circles) 
denote first (second) order phase transitions deduced from QMC 
data. Thin dotted lines are mean-field results. The other lines (solid 
resp. thick dashed for second and first order transitions) are obtained 
by interpolating between numerical data. The data for low densities 
including the paired superfluid (PSF) have been taken from Ref. 13. 



mation (SCA) of a minimal model of hard-core bosons on the 
square lattice defined by the Hamiltonian: 

H = -tY^iMbj+h.c^ +V^mnj - A«E n * 

(i,5) i 

-t'J2 E rn[bi +s b i+s > + h.c] (1) 

i 5=±x;5'=±y 

where m — b\b i £ {0, 1} is the boson density at site i, \i the 
chemical potential, t the n.n. hopping amplitude, t' the am- 
plitude of the correlated hopping, and V the n.n. repulsion. 
The correlated hopping term describes a process where a par- 
ticle can hop along the diagonal of a square plaquette provided 
there is a particle on one of the other two sites of the plaquette. 
In the context of weakly coupled dimers with intra and inter- 
dimers exchange J and J', t can be arbitrarily small, V is of 
order J', t' is of order J' 2 / J, and fx = g^ B H - J + 0{J')[4]. 
Throughout the paper, the energy scale is fixed by t + 1' = 1. 

The case without correlated hopping (t' = 0) has already 
been investigated thoroughly [14]. For strong enough repul- 
sion, an insulating phase with checkerboard (CB) order ap- 
pears at half-filling. The phase diagram is symmetric about 
n = 1/2 in that case due to particle-hole symmetry, and the 
transition from the solid to the superfluid phase is first order 
with a jump in the density [9]. So, there is no supersolid phase 
in the absence of correlated hopping. In the following, we 
study how this picture is modified when correlated hopping is 
introduced. 

Let us first briefly review some technical points. In SSE 
QMC, the various bosonic phases can be determined by study- 
ing the density n, the static structure factor at the wavevector 



(7r, 7r) relevant for the checkerboard solid 

and the superfluid stiffness 

Ps = ^(W* + W*), (3) 

which signals the presence of a SF phase. Here W x and W y 
are the total winding numbers in x and y directions. If both 
order parameters S(ir, ir) /N and p$ are finite in the thermo- 
dynamic limit, the system is in a SS phase. 

The SCA [15] relies on a mapping onto a S = 1/2 spin 
model using the Matsubara-Matsuda [16] representation of 
hard-core bosonic operators = S + , b = S~, and S z = 
1/2 — n . Note that the correlated hopping term transforms 
into a three-spin term. The energy, including zero-point fluc- 
tuations around the possibly non-collinear classical ground 
state [17] calculated within linear spin-wave theory, is then 
minimized assuming a 4-site unit cell to allow for broken 
translational symmetry. The superfluid phase corresponds 
to a ferromagnetic state with a non-zero component in the 
xy-plane, the checkerboard solid phase to Neel order with 
wave- vector (tt, it) and with all spins parallel to the z direc- 
tion, while the supersolid phase discussed below is close to 
Neel order, but the spins acquire a small ferromagnetically or- 
dered non-zero component in the .xy-plane. As a thumb rule, 
the approximation is expected to be qualitatively and semi- 
quantitatively reliable provided the mean-value of the local 
number of Holstein-Primakoff bosons is small compared to 
2S, which is the case here. 

We now explore the phase diagram using both methods 
by first considering a case where correlated hopping domi- 
nates (t 1 = 0.95). The phase diagram in the (1/V, /i/V) 
plane deduced from QMC and from the SCA is summarized 
in Fig. 2, while simulation results for V = 2.2 (1/V = 0.45) 
are shown in the left panel of Fig. 3. Two features of this 
phase diagram are striking: First of all, the large V (small 
1/V) region is dominated by a very large supersolid phase 
that appears for densities above 1/2. Secondly, this super- 
solid phase extends far below the critical value V^ B = 2.38 
(1/V C CB = 0.42) for the development of CB order, down to 
V c ss = 1.7A(1/V C SS = 0.57). In other words, with correlated 
hopping, a supersolid phase can exist without a neighboring 
solid phase as /i is changed. For frustrated quantum antifer- 
romagnets, this implies that supersolid phases can show up 
even in the absence of magnetization plateaus. Note that cor- 
related hopping appears to be crucial for this physics to be re- 
alized. Single particle further-neighbour hopping has been re- 
cently shown to induce supersolid phases near different kinds 
of solid phases [18], but, as far as we can tell, never to stabilize 
a supersolid phase without an adjacent solid one. 

Upon increasing the chemical potential /i, the SF-to-solid 
and SF-to-SS transitions are first order, while the SS-to-SF 
and solid-to-SS transtions are second order, as illustrated in 




FIG. 3: (Color online) Bosonic phases revealed by: a) the density n, 
b) the static structure factor S(tt,tt)/N, and c) the superfluid stiff- 
ness p as a function of the chemical potential for two representative 
cases. Left panel: t' = 0.95 and V = 2.2 (1/V = 0.455). Right 
panel: t' = 0.75 and V = 2.8 (1/V = 0.36) SF stands for super- 
fluid, CB for checkerboard solid, and SS for supersolid. 



FIG. 4: (Color online) Zero-temperature phase diagram for V = 2.8 
as a function of the correlated hopping t' versus the chemical poten- 
tial p. Open squares (closed circles) denote first (second) order phase 
transitions deduced from QMC data. All the lines have been obtained 
by interpolating between numerical data. Thin dotted lines and the 
boundaries to the empty and full system are mean-field results. 



Fig. 3. Interestingly, the SCA and QMC agree on this point, 
as well as on the overall structure of the phase diagram. The 
SCA clearly overestimates the extent of the supersolid phase, 
but it correctly predicts that it extends below Vcb- This is 
an important remark for more realistic models with possibly 
positive correlated hopping amplitudes. In that case, QMC 
will not be possible due to the minus sign problem, but the 
SCA can be expected to be qualitatively reliable. 

Next, we investigate how this large supersolid phase 
evolves from the case without correlated hopping (t' = 0), 
where there is no supersolid at all. We consider an intermedi- 
ate value of the repulsion (V = 2.8) and follow the evolution 
of the T = phase diagram as a function of t'. The results 
in the plane (£', /j.) are plotted in Fig. 4. This phase diagram 
has been obtained with QMC simulations on lattices with up 
to (24 x 24) sites using temperatures (3 = 2L. The finite size 
effects are remarkably small, as illustrated in the right panel 
of Fig. 3, where the results for t' = 0.75 obtained on different 
clusters for n, S(tt, it) and ps are plotted as a function of [i. 

A striking feature of this phase diagram is the strong 
asymmetry introduced by correlated hopping compared to the 
particle-hole symmetric situation when t' = 0. When cor- 
related hopping is introduced, the phase separation above the 
plateau is rapidly replaced by a supersolid phase, which grows 
continuously to extend over all the high density region when 
only correlated hopping is present. In contrast, phase separa- 
tion persists below the solid phase for all values of t' . 

The quantum phase transition between the solid and the SS 
falls into the conventional SF-insulator universality class, as 
expected since the gapped excitations of the solid order are not 
expected to influence the nature of the quantum phase transi- 
tion, and in agreement with the results recently reported for 
a spin model [7]. By contrast, the transition between the SS 
and the SF at zero temperature is first order below t' w 0.25, 



and seems to be continuous above. For the sizes available, the 
extracted critical exponents are consistent with the continuous 
transition being in the 3d Ising universality class. This might 
be a finite size effect though: The long-wavelength gapless 
excitations of the SF phase are expected to change the univer- 
sality class of this quantum phase transition [19], and could 
therefore give rise to a crossover phenomenon at large length 
scales. 

Next we discuss the thermal transitions of the supersolid 
phase. In general one expects two melting transitions for a 
SS: a KT transition when the SF stiffness vanishes, and an an- 
other one when the solid order melts whose universality class 
depends on the type of order [20]. This question has already 
been addressed, first for hard-core bosons on the triangular 
lattice [21], then for a spin model with anisotropic exchange 
integrals. [7] In both cases, two phase transitions have indeed 
been observed. In the case of the spin model, closer to the 
present case since the melting of the solid is in the Ising uni- 
versality class, the KT transition has been found to lie always 
below the Ising transition, suggesting that the supersolid needs 
a solid phase to develop. In the present case, we expect the 
situation to be quite different since at zero temperature a su- 
persolid can exist without a solid. This is confirmed by our 
investigation of the two representative cases for which zero 
temperature data have been shown in Fig. 3: t' = 0.75 and 
V = 2.8 (Fig. 5a and Fig. 5b) and t' = 0.95 and V = 2.2 
(Fig. 5c). For the first case, we indeed find two transition 
lines which both smoothly go to T = upon approaching ei- 
ther the solid phase (KT) or the SF phase (Ising). They cross 
in the middle of the SS phase, defining a region close to the 
solid where the KT transition is below the Ising one, as in 
Ref. 7, but also a region close to the SF where the Ising tran- 
sition is below the KT one. This is even more dramatic for the 
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FIG. 5: (Color online) (a) Thermal melting of the supersolid for 
t = 0.75 and V = 2.8: Binder ratio (empty black symbols) and 
superfluid stiffness (filled blue symbols) as a function of temperature 
for /j, — 6 (upper panel) and fi = 10 (lower panel). The location 
of the melting transitions are marked by vertical dashed lines, (b) 
Finite-temperature phase diagram as a function of the chemical po- 
tential for t' = 0.75 and V = 2.8 (1/V = 0.36). (c) Same as (b) 
for t' = 0.95 and V = 2.2 (1/V = 0.455). In (b) and (c), lines 
are interpolations between numerical data, and error bars are smaller 
than the symbols. 



second case of Fig. 5c, where the melting of the Ising order 
occurs entirely inside the superfluid phase. 

Finally, let us comment on the experimental implications 
of these results. Whenever correlated hopping dominates, one 
can reasonably expect to find supersolid phases, even in re- 
gions where no plateau has been detected. In SrCu2(B03)2, 
such a large domain exists between the 1/8 and 1/4 plateaux, 
in a field range accessible to NMR, a technique well suited to 
detect lattice symmetry breaking. Even if experiments such 
as NMR, which are almost always done in steady field, can- 
not be performed, our results firmly establish the presence of 
two phase transitions, the melting of the solid often taking 
place below the KT transition if correlated hopping is present. 
Since the magnetization is expected to show only a very weak 
anomaly (if any), it would be very interesting to perform sys- 



tematic specific heat measurements of frustrated dimer anti- 
ferromagnets in high field and at very low temperature to try 
to detect ordering in region where the magnetization does not 
exhibit any plateau. It is our hope that the present paper will 
stimulate such investigations. 
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